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Abstract 

We study reflected BSDE in a general continous martingale setting, when the coefRcient 
/ of the driver has at most quadratic growth in the control variable with a bounded 
terminal condition and a lower obstacle which is bounded above. We obtain the basic results 
in this setting : comparison and uniqueness, existence, stability. For the comparison theorem, 
and the special comparison theorem for reflected BSDEs (which allows one to compare the 
increasing processes of two solutions) , we give straightforward, intrinsic proofs which do not 
rely on comparison for standard BSDEs. We obtain existence by using the flxed point theorem 
and then a series of perturbations. While we flrst show existence under the assumption that 
/ is Lipschitz in the primary variable Y ^ we then extend this to the case where / can have 
slightly-superlincar growth and the case where / is monotonous in Y with arbitrary growth. 
We also obtain a local Lipschitz estimate in BMO for the martingale part of the solution. 

1 Introduction 

Since Pardoux and Peng initiated the systematic study of (non-linear) backward stochas- 
tic differential equations (BSDEs thereafter) in (53] , these equations have been proved useful 
for a number of areas. They are intimately linked with the stochastic version of the maximum 
principle for stochastic optimal control problems, they provide probabilistic representations 
for the solution to some partial differential equations (PDEs thereafter), and they are nat- 
ural in finance where, beyond coming in through the use of the above fields, they provide 
the natural language in which the replication of a European option is expressed. Refiected 
BSDEs (RBSDEs) have applications in the same areas, where they are linked to optimal 
stopping, PDEs with obstacle (variational inequalities) and American options. 

"Refiected" BSDEs are BSDEs subject to a constraint : the solution process Y is required 
to remain above a lower obtacle L. In order to achieve this, it is necessary to add to the 
usual dynamics dFs = —fgds + ZgdWs a "force" dK that drives Y upward. One wants 
that extra term to be minimal, so that K is only active to prevent Y from passing below 
the obstacle L. This optimality condition (the Skorohod condition) is therefore expressed as 
/o ^{ys>Ls}^^s ^ 0- in the end a reflected BSDE takes the following form : 

' dY, = -/(s, r„ Zs)ds - dK, + Z,dW„ 

Ft > for alH e [0,T], (1) 

K is continuous, increasing, starts from and / l^Ys>Ls}dI^s — 0, 

Jo 

where the solution to be determined is now the triple (Y, Z, K). In this paper, / will be will 
be assumed to have at most quadratic growth in z. The terminal condition ^ is assumed to 
be bounded and the obstacle process L a continous scmimartingale bounded above. 
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Reflected BSDEs were introduced by El Karoui, Kapoudjian, Pardoux, Peng and Quenez 
in [13]. These authors considered the case where / is Lipschitz, the terminal condition is 
square-intcgrable and the lower obstacle a continuous square-intcgrablc semimartingalc, the 
natural extension of Pardoux and Peng [53] . Given the similarity between reflected BSDEs 
and standard BSDEs. most developments made in the theory of BSDEs can be adapted to 
reflected BSDEs, in particular on the problem of well-posedness. An account of this can be 
found in Bayraktar and Yao [3]. The first existence, uniqueness and comparison results for 
quadratic BSDEs (that is, when the coefficient / is allowed to have a quadratic growth in z) 
were obtained by Kobylanski in |17j . under the assumption that the terminal condition ^ is 
bounded. These results have since then been extented to more general cases, for instance by 
allowing / to have slightly superlinear growth in y in Lepeltier and San Martin [T^] , or ^ to be 
exponentially integrable in Briand and Hu [7] . Those results have been obtained for reflected 
BSDEs as well, see Kobylanski, Lepeltier, Quenez and Torres [T8| and Lepeltier and Xu [2T| . 
Work has also been done regarding the assumptions on the obstacle L, for instance Peng 
and Xu |25| worked with obstacles, and in the context of quadratic / and exponentially 
integrable ^ Bayraktar and Yao [2j removed the condition that L be bounded. However, the 
case of quadratic BSDEs remains signiflcantly more difficult than that of Lipschitz BSDEs, 
and the methods used initially arc often quite involved. Recently, Tcvzadzc gave in |26j a 
new and simpler treatment for the existence of solutions to a quadratic BSDEs, and again, 
more recently, Briand and Elie gave a new approach in [S] . 

Backward stochastic problems in general can be considered outside a Brownian setting. 
For instance, BSDEs have been studied in a general martingale setting in El Karoui and 
Huang [12] , and in a general filtered probability space in Cohen and Elliott [9] . While papers 
dealing with the theory for reflected BSDE are set in a Brownian setting, we consider here 
reflected BSDEs in a general continuous-martingale setting and obtain the basic results 
(comparison and uniqueness, existence, stability) using simple methods, intrinsic to RBSDEs. 

We flrst show, in section [31 how comparison for quadratic reflected BSDEs can be proved 
directly, as in El Karoui et al. [13], using the BMO argument from Hu, Imkeller and Miiller 
[11] . as opposed to via an optimal stopping representation and comparison for BSDEs (see 
[T5]). This result holds naturally for / locally Lipschitz in y, instead of globally Lipschitz. 
We then give a new proof of the special comparison theorem for reflected BSDEs, which 
allows one to compare the increasing processes when one RBSDE solution dominates another. 
This theorem was flrst proved in Hamadene, Lepeltier and Matoussi |10j . and reused in 
Peng and Xu [25. In the papers ([lO], [25], [20], [18]) where it appears, the proof relies 
on the penalization approach and the comparison theorem for standard BSDEs, comparing 
quantities which, at the limit, become the increasing processes. The proof we give here is 
more related to the non-linear Snell envelope view of reflected BSDEs, and holds without 
the regularity assumptions usually made on /. 

In section |¥] we show how the perturbative technique introduced by Tevzadze |26j can be 
adapted to reflected BSDEs to prove existence of a solution. When trying to solve a quadratic 
BSDE with the flxed point theorem, one can conclude only when the data (/(•,0,0),^) are 
sufficiently small (/(s, 0, 0)ds is the residual drift, that drives the solution even if {Yg, Zg) = 
0). One can then built on this result to construct a solution in the general case. 

This technique can be understood as a type of "vertical" splitting and recombination, 
and is in that sense an analogue to what is done for Lipschitz BSDEs. In that classical 
setting, if one works with the natural norm on the space where one looks for solutions, which 
in that context is the space of square-integrable processes, one finds that the fixed point 
theorem applies if the time interval is small enough. A natural way to use this is then to 
split a general time interval into pieces small enough that one can obtain a solution on each 
interval, and patch them together to obtain a solution on the whole interval. For quadratic 
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BSDEs, since one can apparently solve the BSDE only for small data, the idea is to split a 
general set of data into pieces small enough that one can obtain a solution for each piece, and 
then combine them to obtain a solution to the initial problem (see also [IZj). One can also 
understand this method as a series of perturbations. One first solves a BSDE with small data, 
then successively solves pertubation equations and adds the associated solutions, allowing 
the size of the data to grow at each step. At the end, one has built a solution to the initial 
BSDE with macroscopic data. 

In order to obtain a solution to a backward stochastic equation (with small data) via 
the fixed point theorem, one mainly needs to understand the underlying backward stochas- 
tic problem (for BSDEs, this underlying problem is the semimartingale decomposition, for 
reflected BSDEs, it is a Snell envelope problem). However, for the perturbation procedure 
to work well, the underlying problem should be a linear problem (for instance, it has been 
successfully applied recently in Kazi-Tani, Possamai and Zhou [16] to quadratic BSDEs with 
jumps). This way, the equations satisfied by the perturbations are of the same nature as the 
equations satisfied by the solutions. This is not the case for reflected BSDEs. as the under- 
lying problem is an optimization problem, and therefore non-linear. It is however possible to 
identify the equation that a perturbation should satisfy. The obstacle cannot be perturbed 
easily during the procedure, but this can be dealt with by assuming from the start that it is 
negative, a case which covers all the others by a simple translation. In particular, unlike in 
[TS] . we don't need L to be bounded but only require it to be upper bounded. 

We then study the stability of the solution with respect to changes in the terminal condi- 
tion ^ and in the residual drift /(•, 0, 0). This stability can be obtained as a direct consequence 
of the use of the fixed point theorem for showing existence (uniqueness being acquired at 
that stage). We further obtain for the martingale part of the solution a local Lipschitz bound 
in the space BMO. Global Lipschitz bounds in THF were obtained already in Briand, Dclyon, 
Hu, Pardoux and Stoica [S] (see also Briand and Confortola [3], Ankirchner, Imkeller and 
Dos Reis [T]). Kazi-Tani, Possamai and Zhou [12] provide a global ^-Holder estimate in the 
smaller space BMO. Here, we can obtain a stonger regularity for small perturbations. 

In section [5l we extend the scope of the existence theorem of section |4l In the usual 
setting where / is a Lipschitz function of y, the sequence of perturbations described above 
can be performed uniformly without problem. When the first derivative fy is not a bounded 
function of y, the maximal allowed size for a perturbation depends on the size of the solution 
to the reflected BSDE that one wants to perturb, so it is not clear a priori that the procedure 
would terminate after finitely many perturbations. We show, however, that this is the case 
as soon as one can obtain an a priori bound for Y in 5°°, thus recovering the generality 
on / considered in Kobylanski et al. [T^ — slightly superlincar growth — and Xu [55] — 
monotonicity and arbitrary growth. This monotonicity assumption, introduced by Pardoux 
(see for instance [52]) in the framework of BSDEs with Lipschitz-dependence in z, and later 
generalized to coefficients / quadratic in z in Briand, Lepcltier and San Martin [8] , is relevant 
for reaction-diffusion equations. 

In the following section, we motivate the general continuous-martingale setting, and spec- 
ify the notation and the framework that will be used throughout the paper. 

2 Setting 

General, continuous- martingale setting. 

The reflected BSDE ([1]) is set in a Brownian setting. Even in this setting, it is sometimes useful 
to consider that the reference increasing process is not the time ds. For instance, Pardoux and 
Zhang observed in [24] that when looking at the BSDE associated with a scmilinear parabolic 
PDE in a domain with Neumann boundary conditions (so-called "generalized BSDEs"), the 
drift term is of the form f{Ys, Zs)ds + f'{Ys)dAs where A is the local time of the underlying 
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diffusion on tfie regular boundary, and is tlierefore ortliogonal to the Lebesgue measure. In 
that case, one can enhance the increasing process by setting dC ~ ds + dA, and find /" such 
that f{Y„ Zs)ds + f'{Ys)dAs = f"{Y,,Zs)dC (see El Karoui and Huang [H]). 

The reference martingale M doesn't have to be a Brownian motion, and in partic- 
ular doesn't have to enjoy the martingale representation property (see El Karoui, Peng 
and Quenez [2]). The martingale part N of the solution then has the decompotion TV = 
/ ZdM + N-^ on the reference martingal M, with N-^ orthogonal to M (i.e. (iV-"-, M) = 0). 
The quadratic variation d{M) of M is assumed to be absolutely continuous (component- 
wise) with respect to dC (if it is not, one can enhance dC so that it becomes the case). Write 
then d{M) = adC = aa*dC. Since Zg is uniquely determined only when no component of 
d{M)s is zero, the drift term will only depend on Za. The orthogonal component N-^ of the 
martingale part N is meaningful in the context of Follmer-Schweizer strategies in incomplete 
markets (see [H]). One wants in that context to have in the drift a quadratic term gsd{N-^)s. 
By the perturbative nature of the method used later to prove existence, one is naturally led 
to consider also a term d{i', N-^) linear in N-^. 

So in the end, we will study the following general reflected BSDE : 

dV, ^ -dV{Y, N)s - dKs + dNs, 

'--i' (2) 

Yt > Lt for t < T, and ^ ^ 

K increasing, continous, starting from and such that l^Ys>Ls}dKs = 
where the drift is given by 

dViY.N), = f{s,Y,, Z,a,)dC, + d{iy,N^), + g,d{N^),. 
This is referred to as the reflected BSDE of parameters/data (V, ^, L) = (/, i/, g, ^, L). 



The framework is a filtered probability space {^,J^ = (-7^t)te[o.+oo[i P) satisfying the usual 
conditions. T is an J^-stopping time valued in M_|- (finite stopping time). The continuous 
square-integrable martingale M is assumed to be BM O (see below) . All the processes con- 
sidered are continuous. 

C is a continuous and progressively measurable increasing process (starting from 0) 
such that, roughly, all the finite variational processes which are related to parameters (not 
depending on the solution) are absolutely continuous with respect to it. In particular, 
d{M)a = dsdCs = '7a<J*gdCs- It is assumed that the positive symmetric matrix a (or equiva- 
lently a) is bounded away from and infinity (i.e. bounded and uniformly elliptic). 

The parameters of the BSDE (coefficients /, z^, g of the drift V , terminal condition ^, 
obstacle L) arc as follows : 

• / : 17 X [0,r] X R X Mi^d(K) ^ K is Prog ® ^(R) (g) B(Mi.d(M))-mcasurable, where 
M„_d(K) is the space oi n y. d matrices with entries in K. Prog ~ Prog{FT) is the 
progressively measurable sigma-ficld on the interval [|0,r|] (the set of pairs (a;,i) such 
that t < T{io)). 

9 1/ is a BMO martingale orthogonal to M (that is (i/, M) = 0) 

• g is a progressively measurable and bounded scalar process. 

• ^ is an J^T-measurable, bounded random variable. 

• i is a continuous semimartingale bounded above. 

Throughout the paper, we assume that / has at most quadratic growth in the variable z, in 
the following sense : 
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(Aqg) There exists a growth function A(-) (i.e. A : M — > W+ symmetric, increasing on R+, 
bounded below by 1) and a positive process h £ L'^bmo (i-^- / hdM G BMO, see 
below) such that : 

\fit,y,z)\<Xiy){h^,+\z\^). 

The assumption as written above allows for any growth in y, although most of the time / 
will have at most linear growth in y. 

Solutions to the reflected BSDE. 

A solution to the reflected BSDE is generally understood as a triple S ~ {Y, N, K) where Y 
is a semimartingale, N a square-integrable martingale (€ TH?) and K an increasing process 
(starting from 0), such that Q is satisfied, with A = / ZdM + N^. 

Note that a solution can also be understood as a pair S ~ {Y, N) such that, definining 
K from Kq = and the dynamics equation in K is indeed found to be increasing and 
satisfies the Skorohod condition. This will often be what is meant by solution in the rest of 
the paper. 

Under the assumption of quadratic growth and bounded terminal condition, we consider 
only bounded solutions : Y G S°°. For those, N is found to be a BMO martingale (see the a 
priori estimate of proposition [T] below) . So a solution will always be understood as being in 
S°° X BMO ( xA). 

Spaces of processes. 

We will make use of the following spaces. 

• BMO{P) is the space of all the BMO P-martingales, that is those for which the norm 

II^I1Ia/o(p) = sup \\Ep{{N)t - {N)t\moo 

is finite, where is the set of stopping times t such that Q <t <T. The mention of 
the measure P will be omitted whenever no confusion is possible. When X G BMO, 
then £{X) is a UI martingale, so one can define a measure Q by stating that on 
|2 = g{X)t. Also, we will use frequently the fact that for any N € BMO{P), 
N = N - {X,N) is in BMO{Q) (cf Kazamaki [IS], theorems 2.3 and 3.3). 

• Li^BMO ^s space of processes h such that J hdM £ BMO. 

• L°°''^ is the space of processes x such that J^\xs\'^dCs G L°°, and L°°'^ is that of 
processes such that J^\xs\dCs G L°°, with norms ||a;||^_2 = ||/o^|a;spdCs||oo and 
Iklloo.i = ||/o^|a;^|dC,||oo- 

The first space will be used for the process r appearing in many growth assumptions 
on y later on. It is there to take into account the fact that Ct might not be bounded 
(in facts, if C incorporates a local time of a diffusion on a regular boundary, it has 
exponential moments but is not bounded). In the case where dCs = ds and T is a 
constant (or a majorant for the bounded terminal time), one would have r = 1 and 

lkllL,2 = r. 

• ^ is the space of accumulators, that is : progressively measurable, continuous, increasing 
processes starting from 0. 
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BMO property for N. 

Proposition 1. Let f satisfy (Aqg), ly G BMO and g be hounded. Let Y he a continuous 
semimartingale, N he a square-integrahle martingale and K be an increasing process such 
that Y has the decomposition : 

dY = -dViY, N) - dK + dN, 

where dV{Y,N), = /(s, F,, ^sCr,,)dCs + d{u,N^)s + g{s)d{N^)s. IfY is bounded (i.e. Y e 
S°°), then N G BMO and K e Abmo- 



<BMO 



Here, Abmo refers to increasing the processes K & A such that the quantity Iji^'lUf 
supf\\E(^KT — Kt\J-t)\\oo is finite. Note that this statement is, to some extent, not so much 
about solutions to a (possibly reflected) BSDE but about quadratic semimartingales (see 
Barrieu and El Karoui [5]), and quadratic semimartingales are considered here up to a 
monotonous process. 

Remark 2. The result implies in particular the following : if F is a bounded semimartingale 
with decomposition dY = —dV — dK+dN, with K monotonous (which boils down to increas- 
ing, up to considering —Y) and if the process V is in Lgj^jQ (i.e. supji?( |dT4| |-^t) ||oo < 
+oo), then N £ BMO and K e Abmo- 

Proof. The proof uses the usual exponential transform. Let ^ € M, whose sign and value will 
be chosen later. By Ito's formula for the process exp{fiY) between a stopping time t £ 
and T one has 



T 



(3) 



Since Y G S°°, the process e^^ is bounded, and since iV is a square-integrable martingale, 
J e^^ dN is a martingale. We have 

\dV,\ < \f{s,Ys, Z,a,)\dCs + \d{v,N^)s\ + \g{s)\\d{N^) 

and using the quadratic growth assumption on / we have 

\f{s,Y,,Zsa,)\ < \{Ys){hl + \Z,as\^) < K{hl + \Z s<J , 



where A = A(||y||5oo). In order to make the following computations lighter, we use the 

Ft), where ... is replaced 



notation for the weighted squared quantity ■ E(^J^ e^^= 

by the appropriate quantity depending on whether one is measuring h, Z, or and N. 
Using the Kunita-Watanabe inequality and ah < a'^ + b'^, we see that 



T 



E[ I e''^=|dK| J-,) <A + +\^\f. + \N^\f. + \\g\UN 



[t] -r l^l[t] 



|2 I lAr-L|2 I 11 II |Ar_L|2 



RecaU that by the orthogonality of M and TV-^, d{N) = \Za\'^dC + d{N-^), and therefore 
both \Za\'^dC and d{N-^) are less than or equal to d{N). So, setting 6 = (A + 1 + |l.g||oo), 
and taking the conditional expectation of ([3]) with respect to Ft, one has 



- f^\K\lt] + {^- \t^\b}\N\l^ < el^lll^ll + ImI [A\h\ 



0. 
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We now choose fi = —Ah. Since 6 > 1, 6^ > 6. We deal with the weights by saying that 
g-|A'lll"i^ll < gtJ-Ys < glA'llli'll^ Taking the the sup^, we obtain finally 



\K\\A.,.o + \\m 



BMO 



< 



e 



Sb\\Y\\ 



2b 



l + Ab{A\\h\ 



\BMO) 



< +O0. 



□ 



Note however that while it indeed gives a bound for N e BMO, this estimate is not fully 
satisfying in that it does not guarantee that if — > then ||-/V||Bj\fo — > 0. One could 
also provide an a priori estimate for ||K||5oo, but that will not be used in the paper. More 
precisely, the standard assumption under which existence is proved in section U is enough to 
carry the perturbation procedure. In section [5] however, we consider more general cases and 
there we will need the related a priori estimates. 



3 Comparison theorems and uniqueness. 

3.1 Comparison theorem. 

We give below the comparison theorem in our setting, which guarantees uniqueness in 
the existence theorem I 111 The regularity assumption that we require for the theorem to hold 
is, for notational simplicity, the following : 

(Anf ) The function / is of class (in the variable (y, z), for all lu, t) with 

\fyit,y,z)\< p{y)r^ and \Mt,y, z)\ < p' {y)iht + \z\), 

for some growth function p and p' . 

Theorem 3. Consider two sets of parameters (/, i^, ^, i) and (/', i/, 5', L'), and assume 
that : 

1. there exist solutions {Y,N,K) and {Y' , N' , K') to the corresponding reflected BSDEs, 

2. the parameters are ordered : /' < f , g' < g, ^ C '^'^'^ ^ 

3. f is regular enough : it satisfies (Ajjf ). 
Then one has Y' <Y . 

Proof. The proof relies on a classical linearization argument and the properties of solutions 
to a linear BSDE, just like when proving comparison for BSDEs. More precisely, we study 
the positive part (Ay)+, where AX = Jf' — X for a generic quantity X, and show that 
(Ar)+ < 0. 

Denoting by I the local time of ISY in 0, the Ito-Tanaka formula gives 

d(Ar)+ = i{An>o}rfAi; + \dis 

- dAVs - dAKs + dANs 



1 



{AY,>Q} 



i^dls . 
2 
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Now, gathering terms, rewriting differences, and hnearizing some, 



+ d{v, {NY - N^) + 9'sd{iN')^)s - gsdiN"-)., 



dC, 



(A/)(s, n', Z'^as) + fis, n , - fis, n, Z,a,) 

+ d{y, AN^)s + {Ag)s d{{N')^)s + g.s [d{{N')^) - d{N^) 

(A/)(s, y;, + fis, n', - fis, Ys, Z'^a,) + fis, Y„ Z'^a^) - fis, F., 

+ d{u, AN^)s + iAg), d{iN')^)s + gsd{iN')^ + N^, AN^) 



dCs 



Af + FyAY + F.AZa 



dC 



+ d{iy',AN^) + iAg)d{iN')^), 



where 

Fyis) - Fyis, n, y;, = / fyis, + u AY„ Z '^<j du , 

Jo 

F^is) = F^is, Ys,ZsCrs,Z'^as) = / /^(s, Ys,Zsas + uAZs(Ts)du and 

Jo 

u' = u + J gdiN^ + iNY) ■ 
So we can rewrite ^ as 

rf(Ar)+ = -dD - FyiAY)+dC + l{Ay>o} \dAN - F^AZadC - d{v' , AN^) 

where 

dD = l{Ay>o} \AfdCs + Ag d((iV')^> + dAK 
is a decreasing process. Indeed A/ < 0, Ag < 0, > and 

'^{AY>o}dAK = l{AY>o}c^^'' - l{Ay>o}c^^' < 

=0 >o 

because dK > and on {AY > Q}, Y' > Y > L > L' so dK' 0. 

(AF)+ is therefore seen as the solution to a Uncar equation. Define the integrating factor 
Bt = e-^o Py^'^'^'^^" and the measure Q by ^ = £ij F^a'^dM + v')t. By the assumption 
on fz and the fact that h G L^mo ^'^'^ N,N' and are in BMO, j F^cr^^dM + v' is in 
BMO, and therefore Q is indeed well defined. Then, AN = AN ~ J F^AZcrdC, - (i/', AiV-^) 
is a i3MO((5)-martingale. By the assumption on fy, the process -B-1{ay>o} is bounded so 
/ BulAY^>odANu is again a Q-martingalc. Therefore, looking at the dynamic of F = BY 
under Q we finally find that 



-'-dl 



< iAYt)+ = Ec 



< 



□ 
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Remark 4. The previous theorem is stated, for convenience, for a function / which is C^. 
Typically, for the comparison theorem one only requires / to be locally Lipschitz, in which 
case the processes Fy, have to be replaced by the differential quotients : Syf{Y, Y' , Z'a) = 

'^^^ '%-}Zy^'^ ly^y'; etc, and the proof works as long as Fy G L°°'^ and F^ G L^mo- These 
criteria are satisfied as soon as 

(AiLip) There exist growth functions p and p' , and processes r g L°°'^ and h G L^mo such. 
that 

\f{t,y',z') - fit,y,z)\ < piy,y')r^\Ay\+p'iy,y'){ht + \z\ + \z'\)\Az\ . 

Note : when p, p' are constants this is the standard assumption of local Lipchitz regularity 
made in the quadratic BSDE literature (sec for instance Briand and Elie and Tevzadze 
|26j). However, since we are dealing with bounded solutions, the assumption can be weakened 
to the case (AiLip) where p, p' are growth functions. 



3.2 Special comparison theorem. 

When the two sets of parameters arc in a comparison configuration and when the lower 
obstacles are the same, one can say more than Y' < Y and also compare the increasing 
processes of the two solutions, K' and K. 

Proposition 5. Let (/, g, v, ^, L) and (/', g' , v, L) be some parameters, and assume that : 

1. there exist solutions S = {Y, N, K) and S' = {Y' , N' , K') to the corresponding RBSDEs, 

2. the drift coefficients are ordered : /' < f , g' < g, 

3. Y' is dominated by Y :Y' <Y. 

Then it is the case that dKt < dK^. 

The intuition is quite clear. First, since one has F/ < Ft, if y doesn't touch the barrier 
(Ft > Lt), then dKt = and whether F/ > Lt or F/ = Lt, one has dK'^ > = dKt. So the 
only non-trivial case is when F touches the barrier, and therefore F' as well. In that case, 
since the extra forces dK' and dK are minimal, they only prevent the drifts dV' and dV 
from driving the solutions F' and F under the obstacle. But since rfV/ < dVt in that case, 
the correction that could be needed for F will be less than that needed for F'. The proof 
makes this heuristics rigorous. 

Proof. In this proof, contrary to the rest of the paper, AX denotes X — X' for a generic 
quantity X. In order to deal properly with what happens locally when the process AY 
touches 0, we proceed as in El Karoui et al. [13] : write down the structure of AY and AF+, 
argue that these two processes are equal (since by assumption AY > 0), identify their finite 
variational and martingale parts, and then extract the relevant information. Our goal is to 
prove that dAK < 0. 
We have 

dAY = -dAV - dAK + dAN and 
d(AF)+ = l{AY>o}dAY + ^dl , 

where I is the local time of AY at 0. Identifying the finite variational and martingale parts, 
we see that 

-dAV - dAK = l{Ai'>o} ( " dAV ~ dAK^ + ^dl and 
dAN = l[^Y>o}dAN , 
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that is to say 



4Ay=o} 



( - dAV - dAK^ =^dl and 
l[AY=o}dAN = . 



The second equation impUes, by Ito's isometry and the orthogonaUty between M and 
AiV-L, that lAy=o(|AZcr|2dC + d{AN-^)^ = 0. So we know than on the set {Y' = Y} 

(i.e. against l{Ay=o}) we have Y = Y' and Z = Z'. We also notice that by the Kunita- 
Watanabe inequahty, l^^Y=o}d'{'^' , AN-^) = (for any continuous seminiartingale v' that 
could possibly come up). 

The drift term can be rewritten 

dAVt - (/(5) - r{S'))dC + d{v, N^) - d{v, {NY) + 9d{N^) - 9'd{{N')^) 
= (/(5) - fiS') + iAf)iS'))dC + d{v, AN^) + d((0), {N')^) 
'd{N^) d{{N')^)] + {Ag)d{{N')^) 

{f{S) - f{S'))dC + d{v, AN^) + gd{N^ + {N')^,Ah 
\Af ){S')dC + d{{Av), (nY) + (A.9)d((iV')^ 
{fiS)-f{S'))dC + d{iy',AN^)] + [diAV)iS') 

where i^' = v + j gd{N^ + [N')^). By the assumptions on the coefficients, we know that 
d{AV){S')t ~ dit > 0. So we find that against l{AYt=o} we have 



In the end, 



^{AYt=0}dAVt = + l{AYt=0}dIt 



l{AY=o}{-dI-dAK) = -dl , 



so 



'^{AY=o}dAK = "l{Ay=o},i^-2 ^ ^ 



>o 



and so we have proven that l{Ay=o}<^'^^ ^ 0. And when AY > 0, one has Y > Y' > L' — L 
so dK = < dK' , and therefore l{Ay>o}c''^^ !i 0, which completes the proof. □ 



4 Existence. 

We now turn to the existence theorem. We use the pertubation procedure like in Tevzadze 
[26| to construct a solution, under the assumption that the derivatives of / are controlled in 
the following way : 

(Ader) There exists p, p', A > 0, r e L°°'^ and h e L^bmo ^^"^^ ^^^^ 

\fy{t,y,z)\<pr^ and \Ut,y,z)\<p\ht + \z\), 
\fyy{t,y,z)\ < Xr^, \fyz{t,y,z)\< Xrt and \fzz{t,y,z)\<X. 
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4.1 Principle. 

As said in the introduction, the strategy is to first apply the fixed point theorem. To 
perform this, one wiU use only the following assumption on / : 

(AiocLip) The function / is differentiable at (0, 0) (in (y, z), for all (w, s) ), and there exist A > 
such that, writing ps = fy{s, 0, 0) and 7., = fz{s, 0, 0), one has 

- for aU w, 5,2/1,7/2, 21,^2 : 

|/(s, yi, zi) - f{s, 2/2, Z2) - Psiyi - 2/2) - Js{zi - ^2)! 

< >^{rs\yi\ + r^|2/2| + ki| + \z2\j {r^lyi - 2/2I + ki - Z2I) , 

- 7 G ^Ia/o and /3 G L°°^^ (that is : J^lMdCs € L°°), 

which follows naturally from the assumption on the second derivative of / only in (Ader)- 
In all generality, it allows also for quadratic growth in y. So what one actually proves first 
is that when / satisfies this assumption (with possibly quadratric growth in y and z), and 
when the parameters are small enough (in a sense to specify), there exists a solution. 

The perturbations procedure is then carried as follows for a reflected BSDE with obstacle 
L < 0. Split the initial data in n (equal) pieces : for i € {!,..., n}, set — and 
a' = where a — f{-,0,0), so that (^*,a*) is small enough. First, there is a solution 

= (F\ A'l) to the reflected BSDE © with smah parameters (f - a + a^,v, g, ^\ L). 

Now, unless otherwise specified, we denote by x'' the sum E^^j^a;^ , for a general quantity x 

indexed by {1, ... , n}. For i = 2 to n, having obtained a solution S** ^ = {Y^ ^, A^* ^) 

to the reflected BSDE ([2|) with parameters (/ — a + a*"^ , i^, 5, ^ ,L), one incorporates one 
more (a', ^*) in the system. One first obtains the perturbation S"* = (F' , A*, AT') solving the 
pertubation equation 

dY' = -dVHY\N') - dK' + dN' , 

Y^^e , 
^ r^'+Y^>r-'+L^, 

dK'^^ + die > and {dx"^^ + dK'){Y'~^ +Y' > +L')=Q 
with drift given by 

dV\Y\N'),= [fi's'^^ + S')- f{S''^) + a',]dCs+d{v^ J 2gd{W)^,{N')^), + gd{{N' 



Is , 



where i/' ^ = ly + J 2gd{N^ ^)^, and ^ is the translated/recentercd function / around 
S . It satisfies / (0) = so the residual-drift (constant part) in this equation is given by 
a*. So the parameters (/' ^ + a*, g, L) here are small in the required sense. Finally, 
one sums := ^ + S*' to obtain a solution to the reflected BSDE ([2]) of parameters 
{f — a + a\iy,g,^\L). For i — n this provides a solution to the reflected BSDE of interest. 

This allows us to conclude to existence for those reflected BSDEs with negative obtacles. 
But then we can show that up to translation, this covers all the cases where the obstacle is 
upper-bounded. 
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Note that the above perturbation equation ([5|) is not a RBSDE in the variable = 
{Y^ , , K^) because is not necessarily increasing. It could be viewed as a reflected BSDE 
in the variable {Y^,N^,K ) but this point of view will not be used (see the remark after 
proposition [9] and remark [TOl after its proof). Also, note that the solution to the initial, 
small RBSDE can be viewed as a perturbation : S = + 5*^, being the solution to the 
RBSDE of parameters (/ — /(•,0,0),i^, 5,0, L). So it would be enough to study only the 
pertubation equations, but it seemed clearer to treat first the small reficctcd BSDEs and 
then deal with what changes for the perturbation equations. 



4.2 Existence for small reflected BSDEs. 
4.2.1 Underlying problem. 

In order to use the fixed point theorem, we need to check that the underlying problem, 
that is to say the backward stochastic problem that one sees when the drift dVt is a fixed 
process and doesn't depend on the solution, defines indeed a map from S°° x BMO to 
itself. For reflected BSDEs, as was explained in El Karoui et al. [T3|, the solution is the Snell 
enveloppe of a certain process (more precisely, Y+J^ dVs is the Snell enveloppe of dVs). 
We want to check that it lives in the right space. 

Proposition 6. Let V G isA/o (in the sense that sup(|ji5( /j"^|(iV^s||-^t) ||oo < +ooj, ^ G L°° , 
and L be upper bounded. There exist a unique (Y, N, K) G S°° x BMO x A solution to the 
reflected BSDE : 



dY = -dV - dK + dN , 
Y > L and lyyLdK ~ 



(6) 



In particular, this applies when dVs = dV{y, n)s — f{s, ys, Zs<Js)dCs + d{v, n^)s + gsd{n^)s, 
for f satisfying the quadratic growth condition (Aqg), v G BMO, g G and {y,n) G 
5°° X BMO. 



Proof. We know from proposition 5.1 in El Karoui et al. 



that Yt is given by : 



Yt = ess sup E 



dV, + Lrlr<T + £.U=T 



(7) 



where are the stopping times r such that t < t < T, and that {K, N) is the Doob-Meyer 
decomposition of the supermartingale Y + V. Our goal is to check that (Y, N) is indeed in 
S°° X BMO. 

For an upper bound on Yt, we say that 



E 



dVs + Lrlr<T + ^It=T 



for any stopping time r, so Yt < \\V\\ii + ||L+j 



E L 



Tt 



+ l|i+l|oo + 11^+11 

+ ||^^||. For a lower bound, we just say 
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that since Y solves ([51), using the fact that K is increasing, 



T 



Yt=E[^+ / dVs + {KT-Kt) 



T 



>El^^ + J dVs 
so Y is indeed in S°° . 

One can then invoke remark [5] after proposition [T] to conclude that N E BMO. □ 

4.2.2 Existence for RBSDEs with small data. 

First one proves that there is a solution when the data are small and when, essentially, 
the drift is purely quadratic in the solution. 

Proposition 7. LeA A > 0, r > such that r E L°°'^ . Let f satisfy assumption (AiocLipz), 
with parameters (/3 = 0,7 = 0,A, ?■) and be such a ~ /(•,0, 0) G L°°'^ (i.e. : J^\as\dCs S 
L°°). Let V = £ BMO and g be bounded by X. There exists eg — ^o{^,f) > such that if 

E=U\\oo + \\f{;0,0)\\L^.i + \\L+\\^<eo , 

then there exists a solution S = {Y,N,K) € 5°° x BMO{P) x A to the reflected BSDE (0j 
with parameters (F, ^, L), where dV{Y, N)^ = f{s, Yg, ZsCrs)dCs + gsd{N-^)s . 
More precisely, 

eo(A,r' 



2iOA(|| 



2 



Also, for any R < i?o(A,r) = — — / ^ ^ — y — ^^^''^ ^^^'^ solution is known to 



satisfy 



Proof. We study the map Sol : S°" x BM O — > S°° x BM O which sends (y, n) on the solution 
(y, TV) to the reflected BSDE 

dY ^ -dV{y, n) - dK + dN , 

Yt^C, (8) 
r > and l{Y>o}dK = , 

where dV{y,n)s = f{s,ys,Zs<Js)dCs + gsd{n^)s . This map is well defined according to 
proposition [HI and (F, N) G x BMO is a solution of ^ if and only if it is a fixed point 
of Sol. It will be seen that Sol is not a contraction on the whole space, but it is on a small 
ball, and it stabilizes such a small ball if the data are small enough. Therefore there exists 
at least one fixed point in the space. 

We study first the regularity of Sol. Take s = (y, n) and s' = (y', n') in S°° x BMO, write 
S = Sol{s), S' = Sol(s'), and Ax = x' — x for a generic quantity x. The semimartingale 
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decomposition of AY is dAY = ~dAV — dAK + dAN , and the terminal value is 0. There- 
fore, applying Ito's formula to {AYY between t S and T, and taking the expectation 
conditional to Ft one has, since AYdAN is a martingale, 

{AYtf +e(^J d{AN)s\J't) =0'- + 2e(^J AY,dAV,\J^t) (9) 

+ 2i;(^y^ AYsdAK.lTt'j - . 

Let us now look at the third term on the right-hand side. Using the fact that YdK = LdK 
and Y'dK' = LdK' one has 

AYdAK = (y' - Y)dK' - {Y' - Y)dK 

= {L- Y)dK^- {Y' - L)dK^ < . 

Let us now deal with the second term : 



e(^J^ AY,dAVs\Tt) <\\AY\\^e(^J^ \dAVs\\Tt 



The assumption on / gives 



\dAV,\ < \[r,\y,\ + Tsly'J + |z,,cr,| + \z[^a.,\j {rs\Ay,\ + \Az,a,\)dC, 
+ \gs\\d{An,n + n')s\. 

Using the Cauchy-Schwartz and the Kunita-Watanabe inequalities, and the elementary in- 
equality (X]r=i ^i)'^ — "■X] have 



E(^j^\dAV\ 



< 2ix{\\rf\\yT + \\rf\\yf + pJaf + \\zaf) ' {\\rf\\Ayf + \\Azf) ' 
+ l|.g||||An^||l|(ny . 

By orthogonality, one has ||z|| and smaller than ||?7.||. By definition of the norm on 

S°° X BMO, + ^ ||s||2. This in particular implies that ||n|| < ||s||. So, using 

< one finally has the estimate 

i?(^y|''^MAy||j-,) <2iA(||r|p + l + l)(||.'|| + ||.l|)||A.|| . 

Equation ^ then yields, using 2ab < + 4&^, 

(Ay,)' +e(^{AN)t - {AN)t\Tt) < \\\AYf 

+ Ax2^X^{\\rf + 2fx2{\\s'f + \\sf)\\Asf , 
and by taking the sup, we finally find that, since ||Ay|| < HAS"!!, 

||A5|p <28A2(||r||2+2)'(i|.s||2 + j|.s'!|2)||Asjl2 . (10) 
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Let's now study the size of S* = Sol(s). Following the very same computations and argu- 
ments as for AS we have first 



(11) 



Since YdK = LdK and, importantly, since K is increasing, one can write 



YdK 



LdK < \\L+\\ Kt -K, 



dV + {Nt - Nt) 



so that 



El^j\dK\F^ < \\L+\\\Yt\ + \\L+\\U\\ + \\L+\\El^£\dV\ 
Reinjecting this into (fTTj) leads to 

iYt)^ + E(^{N)tMTt) < (2l|^f + 10||L+||2) +l||rf +9i?(^ 



\dV\ 



Now, by the assumption on /, 



\dV\< f{s,0,0) + X{rs\ys\ + \zsas\y dCs + \g\d{n) 



so 



E 



One obtains in the end 



\dV\ ^t] < ||/(.,0,0)|U,i + 2A(||ri|2+2)||s|p 



<29£;2 + 2SA2(||r||2 + 2)'i|.sf , 



(12) 



where i?=||e|| + |li+ll + ||/(-,0,0)||oo,i- 

To have Sol be a contraction on a closed (and therefore complete) ball B{0, R) of 5°° x 

BMO, we see from (HUl) and dH]) that we would like R and E to be so small that 2^A^(||r||2-|- 
2)^i?2 < i(< 1) and 2^E'^ + 2^X'^{\\r\\'^ + 2)^i?^ < i?^. This is the case as soon as 

R < i?o(A,r) 



25A(||r||2 + 2) 
i? Ro{X,r) 



□ 



Wc now remove the assumption that the linear terms in the drift are null. 

Proposition 8. Let X > 0, r > such that r E L°°'^ . Let f satisfy assumption (AiocLipz); 

with parameters {l3,j,X,r) and be such that a = /(•,0, 0) G L°°'^ (i.e. : J^\f{s,0,0)\dCs € 
L°°). Let V e BMO and g be bounded by X. There exists eg ~ £o(/3i A,r) > such that if 

E^U\\oo + \\f{;0,0)\\L^.i + \\L+\\^<eo , 

then there exists a solution S = (Y,N,K) e S°° x BMO{P) x A to the reflected BSDE (0j 
with parameters (/, g, ^, L) . 
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More precisely, 

Also, for any R < Rr)(X,r) = — ttt—^ where X = eJI'^HA, if E < ^—l^^, then this 

'J y _ m , y 25A(||r||2+2) ' - ^ ' 

solution is known to satisfy 

ll^llQ = ril|.^ + ll^!lL/0(Q) , 

where Yt = e^o PudCuy^ and N is the martingale part ofY under Q : ^ = £(J ja^^dM + 1/). 

Proof. Write f{t, y, z) = Pty + 7*2: + h{t, y, z), where Pt = fy{t, 0, 0) and 74 = f^it, 0, 0) (so 
that h{t, 0, 0) = /(i, 0, 0) = at). Note that h satisfies (AiocLipz) with parameters (/3 = 0, 7 = 
0,A,r). 

The idea is that if {Y,N,K) is a solution to the reflected BSDE ([2]), one can eliminate 
the linear terms {PtYt + 'ytZt<Jt)dCt + d{v,N^)t in dV{Y,N)t by a pair of transforms and 
obtain a reflected BSDE with purely quadratic drift. Proposition [7] guarantees the existence 
of a solution to such a RBSDE. so undoing the transforms yields a solution to 

More precisely, define the measure Q by |p = £{L) = £{J ja~^dM + v). Then M := 

M - (L, M) ^ M - J ja*dC is a SMO((5)-martingale. Define also B = cxp PudCi^ , 
which is a bounded process. Define 

Hs,y,z) = Bsh{s,B~^y,B^^z) , 
9s = B~'^9s , 

L = BL . 

Note that /i satisfies (AiocLipz) with parameters (/3 = 0,7 = 0, A,r) where A = Ae'l'^IL 

Proposition [7] ensures the existence of a solution {Y,N, K) e 5°° x BMO{Q) x ^ under Q 
to the reflected BSDE © with parameters [h, ly = 0,g, ^, L). 
Indeed, 

llfflloo < exp(|l/?|l)|lg|l < exp(|l/3|l)A = A < +00 , 
||/J(-,0,0)|U,i <exp(||/3||)||/(-, 0,0)11 <+^, 
ll^loo <cxp(||/3||)!|e|| , 
||i+||co <cxp(||/3||)||L+l| < +00 , 

so if _E < exp(— ||/3||)eo(A, r) = exp(— 2||/?||)eo(A, r), proposition [7] applies. 

Define then F = B-^Y, N = gB-^dN = / ZdM + and A' = J^B-^dK. The 
Girsanov (Q — > P)-transform of A^, 

iv = TV - (Z,iv) 

= J ZdM + N^+ J -fZadC+{v, N^) 
= J ZdM + , 
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is a i3MO(P)-martingale. F is a bounded semimartingale, since B ^ is bounded, and differ- 
entiating Y = B^^Y shows that {Y,N,K) is solution to the reflected BSDE Q with data 
(under P). □ 

4.3 Perturbation of a reflected BSDE. 

We now deal with existence for perturbation equations like ([S]). We assume we have a 
solution S*^ = (y^,A^^,A'^) to a reflected BSDE, and want to construct a solution S to 
a reflected BSDE with slightly different parameters. The idea is to construct the difference 
S'^ = (Y"^ , N"^ , K^) = S — S^. The next proposition shows how this can be done despite the 
fact that is not an increasing process, so long as one docs not change the obstacle. 

Proposition 9. Let f satisfy (Ader) with parameters (p, p' , X,r, h) and be such that a = 
/(•,0,0) G let g e be bounded by \ and v ^ BMO. Let also G L°° and be 

upper bounded. Assume that there exists a solution 5^ = (Y'^,N'^,K^) to the RBSDE (0j 
with parameters {f,g,v,£},L^). Now let G L°° and o? G L°°-^ (and = 0). If 

e-2plkf 

then there exist S'^ = (F^A^^ A'^) where G , G BMO{P) and A'^ has finite 
variation, solving the perturbation equation 



' dY^ = -dV^{Y^,N^) ~ dK^ + dN^ 

ys = e , 

Y^ +Y^ > L^ +L^ , 



(13) 



with drift 



dK^ + dK"^ > and l{Y^+Y2>L^+L^}{dK'^ + dK'^) = 



dV^ = [/(s, + Y}, Zla, + Zla,) - f{s, Y^^Z^a,) + aj] dCs 
+d(v+ [ 2gd{NY,{NY)^+9sd{{NY)^- 



2 

So S := S^ + S"^ is a solution to the RBSDE (0j with parameters [f + a^, g, ^^ + ^'^, L^) 
We further la 

solution satisfies 



We further know that for any R < Rn(2X,r) ~ — , ^ r-, ifSE < — R, then this 



\\s%^\\Ym^ + \\Nm,,oiQ)<R' ■ 

Note that while S'^ = (F^, N^,K^) is not the solution to a reflected BSDE, {Y^,N^,K^) 
is. However the drift there would be dV^ — dK^, whose residual action (when (F^, N^) ~ 0) 
is a^dCs — dKl, and this has no reason to be small. We can therefore not simply invoke 

2 

proposition |S] for (F^, A^^, AT ) and have do some extra work. 

Proof. The majority of computations that would need to be done here, related to the dy- 
namics of F^, are very similar to those in the proposition [5] about small RBSDEs, so we only 
do the part which is different. 

Define /(s, y, z) = f{s,y + Y},z + Zl<7s) — f{s, Y^ , Z^ag) + a'^. Note that since / satisfies 
(Ader), / satisfies (AiocLip) with parameters {/3,j,2X,r), where /? = fy{-,Y^ , Z^a) and 
7 = /,(-,Fi,ZV). We have ||^||oc,i < p\\rW' < +oo and 7 G L|mo- 
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Following the same approach as for RBSDEs, we first look at the underlying perturbation 
problem of finding {Y^jN^^K"^) solving the perturbation equation (jl3[) when 

+d(i. + J 2gdiNY,{nY)^+9sd{inY)^ ■ 

Because = + is then solution to the reflected BSDE ^ (or rather ©) with drift 

dvl = dVHY\N^)s + dVHy\o^), = [f(s,y^,+Y,\ z^,as + Zla,)+a^,]dCs+d{iy, +n^),+ 



sd{N^ 



^2 



' n 



5, proposition [5] guarantees the existence and uniqueness of S and therefore of 



S . This allows to define a map Sol' from S°° x BMO to itself. 

Now, to find a solution to the perturbation equation (|13p . we proceed like in propo- 
sitions [7] and |51 the difference being in dealing with dK^ which is not monotonous anymore 
here. Up to doing the usual transformations (proposition [8]) , let us assume that the drift is 
purely quadratic as in proposition [T] Then, Ito's formula first leads to the estimates 



\AYt^\^ + E{{AN^)t,T 



Tt < 2E 



\Y,'\' + E{{N')t 



2E 



Y^dV? 



2E 



2E 



AY^dAKl 

T 

Y^dKt 



For the term in AY^dAK^ one has (even if ^ 0) 
AY'^dAK^ 



(y2' _ 


-Y'^)dK^ - 


- (y2' _ 


- Y^)dK^ 




O 9' 

Y ){dK 


-dK^) 


- r ' - 


r'- 


-Y )dK - 


- r ' - 


9 O 

- Y )dK 




9 9' 

Y )dJ<^- 




9 9 

L)dK^ 


— ^ 

<o 


>o 


>0 


>o 



dK^) 



For the term in Y^dK^ one has however, since 
Y'^dK'^ = (y2 - L2)dii'2 + L^dK'^ 

L^))dK 



0, 



((F' - yi) - (I' - L^))dK^ + L2dif2 



= ((y' - l') - (yi - L^)){dK^ - dK^) + L^dK^ 

= (y' - T'jdK^ - (y' - L')dA-i - (yi - L^)dK'' 



(yi - L^)dK^ +L^dK^ 



>0 



>0 



< L^dA'^ = 0. 



Having observed this, the rest is like the analysis of the map Sol and the eq is the same. 



So in the end, provided that that 5E = ||^^||oo + ||q;^||l=c.i < eo(p, 2A,r) 
there exists a solution (y^, A^^, A"^) to the perturbation equation (|T3)) . 



2iO(2A)(|lrP+2) ' 
□ 

Remark 10. Note for later use that uniqueness holds for the perturbation equations. First, 
under (Ader), (AiLip) holds and so does uniqueness for reflected BSDEs. Then, one can 



—2 



y2 and y'^ 



yi + y 



'2 



argue that if Y^ and Y ^ are two solutions to (fT3)) . then Y = Y 

2 ' 2 / 

are two solutions to the same reflected BSDE, so y =Y and therefore y^ = y ^. Or we 

2 

can also argue that if (y^, N'^,K'^) is a solution to ([T^ . then (y^, N'^,K ) solve a reflected 
BSDE ([2]) for which uniqueness holds. 
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4.4 Existence theorem. 

We can now prove the existence theorem. 

Theorem 11. Let f satisfy (Ader) with parameters (p, p', A, r, h) and be such that /(•, 0, 0) £ 
L°°'^. Let V G BMO, g £ L°° be bounded by \, £^ £ L°° , and L be upper bounded. There exists 
a solution {Y,N,K) £ S°° x BMO x A to the RBSDE (0j with parameters (/, g, i/, ^, L). 

Proof. The proof is done in two steps. First, we show that one can indeed reduce the problem 
to the case L < 0, by translation. Existence for the RBSDE with L < is then proved by 
repeatedly perturbing a solution to a similar RBSDE with smaller data. 

Step 1. If {Y, N, K) is a solution to the RBSDE, and U is an upper bound for L, set 
1^ = y - [/. We see that 

' d~^ = dY ~^dU = -dV -dK + dN -0 = -rfT^ - dK + dN , 

Y = Y -U >L-U =: L , 
. ^{^>L-u/^ = ^{Y>L}dK = . 

Here 

d^} = dV{Y, N) 

= dV(^ + U, N) 

= /(s, ?" + [/, Z<j)dCs + d{v, N^) + gd{N^) 
= Za)dCs + d{iy, N^) + gd{N^) . 

It is clear that ~^ still satisfies (Ader) with parameters (p, p' , A, r, h). And from the assump- 
tion on fy one has 

0,0)1 = [7,0)1 < \fis,0,0)\+prlU , 

so 'a' = 7'(-,o,o) £ L°°'\ 

In the end, (Y , N, K) £ S°° x BM O x ^ is a solution to the reflected BSDE of parameters 
{ J g, ^ , L) satisfying the same assumptions, but with L < 0. 

Step 2. We now focus on the case L < 0. Consider eq given by proposition 1^] (or proposition 
[H it's the same). Choose an integer n big enough so that, setting (where a = 

/(•,0,0)) and e = one has = ^ + \W% = iE<e^. 

First, by proposition [51 there exists a solution (y\iV\i^i) to the RBSDE Q with 
parameters (/ — a + a^, g, i). 

Next, for i = 2 to n, having obtained a solution (F'"\lV*"\^"^) to the RBSDE 

with parameters (/ — a + a'"^, i', g, ^ , i), proposition [2] provides a solution (y% A^*, if') to 

the perturbation equation ([5]) and therefore a solution (F*, A^*,if') to the RBSDE ^ with 

parameters (/ — a + a*, i^, ^ , L). For i ~ n, since ^ = ^ and c?'^ = a, (F , Z , K ) is a 
solution to the RBSDE of interest, which ends the proof. □ 
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4.5 Stability in 5°° x BMO. 

Given that uniqueness holds, the a posteriori bounds that come with the construction of 
a perturbation S" — S' to a solution S in proposition [H] readily give the continuity of the map 
(C, a) ^ (F, N), from x L°°'i to 5°° x BMO. 

We now derive an estimate which shows that it is locally Lipchitz, by a sort of bootstrap 
argument on the above stability result. We use the notation 5S' ~ S' — S, SS" — S" — S for 
perturbations around a fixed solution S and AS* = S" — S" = 6S" — SS'. 

Proposition 12. Suppose that f satisfies (Ader) with parameters (p, p' , X,r,h), that a = 
/(■,0,0) e L°°'^, that V G BMO, that g is bounded by A and that L is upper bounded. We 
consider ^ S L°° and the solution {Y, N, K) to the reflected BSDE of parameters (/, v, g, L). 

Now, for any and {i,",5a") e x let S' = {Y',N',K') and S" = 

(Y" , N" , K") be the solutions to the reflected BSD Es of parameters [f + 5a' ,v, g,£f , L) and 
{f + 5a",v,g,C',L). If 

5E',5E" = \\5^-\\ + \\5a\\<^eo{p,2\r) ^ 



V2 ' ' V2 2io(2A)(|lr||2+2) ' 

where /? = fy{-,Y, Za), then we have 

- < 25e2|l^ll die" - C'lloo + \W' - a'||oo,i) and 

\\N" - N'\\bmo(P) < 2'e'M\c{Y, N){U" - C'llco + - , 
where C{Y,N) is a constant depending on {Y,N). 

Proof. We know that /(s, 5y, 5z) = f{s, Ys + 6y, ZsO-g+Sz) — /(s, Kj, ZgO-s) satisfies (AiocLip) 
with parameters (/3,7, 2A,r), where /3 = fy{-,Y,Za) and 7 fz{-,Y,Za) satisfy ||/?|ioo,i < 
p||r|p and 7 £ L'^bmo ' ^^'^ ^ ^ ^ + / '^gdN £ BMO. We linearize / like in proposition [5] : 
7(s, 5y, 5z) = PJy + jjz + h{s, 5y, Sz). 
Since AY — Y" — Y' has the dynamics 

dAY, - - [Aa, + '^^AY, + 7,,AZ,a, + (his, 6Y^',6Z'Jas) - his, SY' , SZ'^a,)}] dC, 
- d{V, {AN)^)s - gsd{{5N")^ + {5NY, {AN)^)s - dAK, + dAN, , 

doing the usual transformations, with ^ = £{J ^a^^dM +77) and B = e-'o and the 

standard computations give, like in pop . 



|A5||^<2«A^ + 292A (WrW' + 2)' {\\dS^% + \\5S%) .\\AS\\l 



But SS" and 5S' arc the unique solutions to the perturbation equations with parameters 
(/ + 6a" ,T',g,S£_" ,L) and (/ + da' ,V, g,S^' , L), and by the way they were constructed in 

proposition ini (recall that 5E' , 5E" < ^-55^-^^^^^=^) we know that they satisfy 

ll^'^'ll^, \\SS%<^^^^^, so 

II^'IIq + < i?o(2A, r)2 = ^ . 

2102A (||ri|2 + 2)2 

Reinjecting this in the previous estimate we have || AS'||g < 2'^A£' + |||AS'||q and therefore 

|A51|^ = IIArill^ + ||AiV||^,,o(Q) < 2'°A^ 
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Then this imphes that ||Ay|| < 2^ AS and so ||Ar|| < 2^621"'" A£'. For the same reason, 

||AiV||BMO(Q) < 25e2irai AS. By theorem 3.6 in Kazamaki, || AiV||BMO(P) < C{Q)\\Kn bmo(Q)\\ 
where the constant depends only on Q, or equivalently on the martingale J ja^^dAI +V, 
and in fine on {Y,N). □ 

Note that the interesting part of the above result is only the martingale estimate. Indeed, 
the estimate for Y" — Y' in S°° actually holds for any size of parameters (as can be seen easily 
by linearizing the drift, doing a change of measure to get rid of all the terms in N and solving 
for y). As mcntionncd in the introduction, we know that (^, a) i— )• iV is global Lipschitz in 
Ti^, and ^-Holder in BMO. The above estimate shows it is in fact locally Lipschitz in BMO. 

5 Existence under more general assumptions. 

In sectional the existence of a solution was proved under the usual assumption that / is 
a Lipschitz function of j/, and therefore at most linear in y, while being at most quadratic in 
z. In this section, we deal with more general assumptions on /. 

To some extent, we would like to replace the constants p,p',X in (Ader) by arbitrary 
growth functions (while, of course, still assuming that / ends up with a growth in y com- 
patible with existence of solutions). Looking back at proposition 1^1 we see that when p is a 
growth function, the maximal size e allowed for a perturbation (S^, Sa) of the parameters 
would depend on the size of the solution. It is therefore not guaranteed that one can 

choose eo uniformly for the perturbation procedure in the proof of theorem or to put 
things differently, that a sequence of pertubations could terminate in finitely many steps. 
This however can be guaranteed if one can obtain an a priori bound for the solutions Y^s of 
reflected BSDEs with drift {f,iy,g). 

Superlinear growth in y. 

We show in the following proposition how this would work when / is allowed to have slightly 
superlinear growth in y. 

Proposition 13. Consider a set of parameters {f^v,g,S,,L) satisfying the assumptions of 
theorem \lll but with p' , X in (Ader) being growth function instead of constants. Further as- 
sume that\f{t,y,0) \ < \f{t^Q,Q)\+ip{y) for a growth function if such that °° ^^^dy ~ +oo. 
Then there exists a solution {Y, N, K) to the reflected BSDE with parameters (/, v, g, ^, L). 

Proof. We will apply the perturbation procedure as was done previously when p, p' , X were 
constants. 

First, by the estimate in theorem 1 in Kobylanski et al. [18] . we know that there exists a 
function F increasing (a growth function) such that for any set of parameters {f,v,g,^,L) 
satisfying the assumptions and for any solution {Y, N, K) we have ||y|| < S(||^|j, ||a|j). Now, 
for a fixed set of parameters, we define pm — p{F{\\£,\\, We fix n big enough that 

^ p-2p^\\rf 

n - 2io(2A(l))(||rP + 2) = '^^'^'^^ ' 

We will construct n solutions 5* of reflected BSDEs or perturbation equations such that for 
each equation, the size of the data is ^ and the size of the solution is such that < 1. 

Note that the ^ here indicates the multiplication by ~ exp(/g fy{S'^u ^)dCu). 

We know that we can do a transation to be reduced to the case L < so wc assume from 
now on that L < 0. Define C = -f and a* = -a. 
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For i = 1, we first build a solution = (Y'^,N'^,K^) to the reflected BSDE ^ with 
parameters {f — a + , ly, g , , L) . Proposition |8] as it is stated doesn't strictly apply, but we 
can adapt the proof. We define the integrating factor B = ^'^'-^ with /? = /? = (•, 0, 0) € 
L°°'^ and the new measure Q by ^ = £{J -fcj-^dM + v) where 7 = 7° = fz{-,0, 0) G L|^/o. 

Then, like in proposition [71 we look for a solution {Y^,N^,K^) to the reflected BSDE with 
no linear term by the flxed point theorem. We look for a solution in a ball of radius R and 
now further demand that i? < 1, so that the conditions to be met are that 

R < Ro&m^r) = ^ and E, = ^ < ^"^^^^^'"^ = eo(2Aa),0 , 

252A(l)(||r||2 + 2) " 25 

where 2A(1) = e"^ ll(2A(l)). Now, since we have chosen n such that — < — r-, 

y ' y ^ " ^ n - 2iO(2A(l))(||r|P+2)' 

" " " 2io(2A(l))(||r||2 + 2) " 2i02A(r)(||rP + 2) " 2i02A(r) (||r||2 + 2) ^ 

because ||/3 || < p(0)||r|l2 and p{0) < p„. So we indeed get a solution (V^, N^, K^) and doing 
the reverse transforms gives a solution = {Y^ , , K-^) to the reflected BSDE with the 
linear terms. 

2 — 1 

For i = 2 . . . n, we have a solution 5* to the reflected BSDE ([2]) with parameters (/ — 
a + c?'~^, v,g,^ , L) and want to construct the appropriate perturbation 5*. We simply do 
the same computations as in proposition |9l using the integrating factor = ^ '^'-^ 

where = fy{-,Y^ , a), and the change of measure ^ = £{J ^^^a^^dC + V^^^) 

where 7^-1 = /^(•, y'"\ Z'^V) and F^-^ = v + j 2gd(N'~^)^ . 

Because we know by the a priori estimate on solutions of the reflected BSDE that 
lir-'ll < F{\\^-^i\\y^a\\) < nilClUI«ll),wc know that H^"'!! < p(l|F^"'ll)lk|P < 
/5,„||r|p. Therefore, just as above, we find that the size Ei of the data is indeed small enough, 
and so we can produce a perturbation S"*. □ 

As can be seen from the proof, the key to the generalization is to have an a priori estimate 
||y|j5oo < F(||^||oo, ||a||oo,i) foi" some growth function F. 



Monotonicity condition with arbitraty growth in y. 

We can also generalize the result of theorem [TT] to the case where / is so-called monotonous 
(or 1-sided Lipschitz) in y, with arbitrary growth. 

Proposition 14. Consider a set of parameters {f,i',g,^,L) satisfying the assumptions of 
theorem ] 111 but with p,p',X in (Ader) being growth function instead of constants. 

Further assume that |/(t, y,0)| < \f{t,0,0)\+ip{y) for a growth function ip and that there 
exists a constant n such that for all y, y', z, s, uj, 

{y' - y){f{s, y', z) ~ /(s, y, z)) < fir^y - y\^ 

Then there exists a solution (Y, N, K) to the reflected BSDE (0j with parameters (/, z^, 5, ^, L). 

As seen above, it is enough to have an a priori estimate for Y . One can invoke the one 
obtained in the proof of theorem in 3.1 in [28j . Alternatively, having argued that it is enough 
to study the case where the obstacle is negative, one can linearize the driver in the N variable, 
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and do a measure change. Then, using Ito with | P to take advantage of the monotonicity 
condition one could conclude that 

< 2e«^(||e||L + 2||a||Li) = miClMklD' 

6 Conclusion. 

We studied quadratic reflected BSDEs in a general continuous-martingale setting, and 
obtained the basic theorems (comparison and uniqueness, existence, stability) in that frame- 
work. 

We showed how the comparison theorem was easily obtained via the standard linearization 
technique for BSDEs, and under only a local Lipschitz assumption on /. We also gave a 
new proof to the special comparison theorem for reflected BSDEs, which does not rely on 
comparison for BSDEs (through the penalization approach) and therefore does not require 
smoothness of /. This is understandable as this theorem is rather linked to the minimality 
condition of nonlinear Snell envelopes than on a regularity property. 

We then adapted to reflected BSDE the technique introduced by Tevzadze in [5^. The 
first part, where one uses the fixed point theorem to construct a solution for RBSDE with 
small data, essentially controls the superlincarities by imposing a priori some bounds on 
the sought solution. It would work for any kind of supcrlinear growth of / in the variable 
y. It would even hold for superlinear-in-z coefficients / if fz{-,0,0) G L^i\jo in fsiCt. The 
perturbation part required a special attention since perturbations to a rcfiected BSDE don't 
satisfy a reflected BSDE (this is linked to the underlying problem for reflected BSDEs, a Snell 
envelope problem, being non-linear ; the technique would work more easily for non-linear 
equations for which the underlying problem is linear). The problems with the difference of 
increasing processes not being increasing can be avoided if one does not perturb the obstacle. 
We therefore applied first a transformation to the RBSDE to be led to study only the case 
where L < 0. We expect that this approach would also work for doubly reflected BSDEs if 
they can be reduced to studying the case where the lower obstacle L is negative and the 
upper obstacle U is positive. 

We then obtained a local Lipschitz regularity estimate for the martingale part N of the 
solution, in the space BMO, which improves on the previously known regularity. The key 
in obtaining this estimate is essentially a bootstraping of an existing stability result for 
(Y, N) G S°° X BMO. While we used in proposition [T^] the a posteriori estimate coming 
from using the fixed point theorem, the proof could be adapted to make use instead of an 
a priori estimate such as the one given in Kazi-Tani, Possamai and Zhou |16| . or even just 
the qualitative continuity in S°° x BMO. 

When / is Lipschitz in y, the series of perturbations can be carried uniformly without 
using an a priori bound on Y , since it was assumed that fy was bounded. We showed that 
the perturbation procedure could still be used to construct solutions when / is not Lipchitz, 
so long as one can get an a priori estimate for |li^||oo- This allowed us to recover the cases 
where / can be slightly supcrlinear in y or monotonous with arbitrary growth. 

Finally we note that the essential ideas here are not specific to backward stochastic 
problems, and could apply to a certain number of nonlinear equations. 
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